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•  Objectives 

•  Jet  in  Crossflow  characteristics 

9  ONERA  Air-to-air  Experimental  configuration  and 
Computational  Setup 

9  AFRL’s  water-to-water  Experimental  configuration  and 
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9  Investigate  mixing  characteristics  of  propellants  in  combustors 
•  CFD  support  for  the  in-house  Themis  project 
9  Provide  numerical  data  for  theoretical  developments 
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Project  Themis 


*  In-house  combustion  device  research  within  the 
Combustion  Devices  Branch  (RQRC)  of  AFRL 

*  Focuses  on  investigation  of  LOX/Hydrocarbon  high 

*  pressure  combustion  devices  through: 

*  Theory  development 

*  Modeling  and  simulation 

*  Subscale  experimentation  at  reacting  and  inert  conditions: 

•  Design  &  test  a  10k  Ibf  combustor  for  liquid  propellants 

•  Requires  good  mixing  of  variable  density  propellants 

-  •  Temperature  uniformity 

•  Concentration  uniformity 


Jet  in 
Crossflow 
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Staged  Combustion  Cycle  Rocket  Engine 


Fuel  Pump 


Turbine 


Oxidizer  Pump 


Nozzle 

Staged  Combustion  cycle 


Diluent 
Injections 


Turbine 


Mixer 


Temperature  Distribution 
On  the  Turbine  Blades 
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Jet  in  Cross  Flow  Characteristics 


Complex  flow  with  multiple  engineering 
applications: 

•  Gas  turbine  combustors 

•  Film  cooling  of  turbine  blades 

•  Air  pollution 

•  Chimney  emissions 

•  Industrial  burners 

•  Chemical  mixing 

•  Wastewater  discharges 
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.  r  Jet  in  Cross  Flow  Characteristics 

VP 


Counter-rotating 
vortex  pair 
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Conduct  a  validation  study  for  a  baseline  case: 

•  Experimental  data  available 

•  Unity  density  ratio 

•  Multiple  confined  transverse  jets 

•  Single  phase/component 


ONERA  experimental/LES  studies  of  an  eight  jet  mixing 


chamber 


Secondary  flow 


Glass  pipe  Plenum  chamber  Jets  flow  entrance 

for  the  jets  flow  In  the  plenum  chamber 
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V  Previous  Numerical  Work  at  AFRL: 

Air-to-Air  Experimental  Configuration 

Diluents  Flow 


i 
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A8r4o*A8r  Numerical  and  Experimental  Comparisons 

Diluents  Flow 


Axial  velocity  Contours,  x  =  2d 


Axial  velocity  contours,  x  =  5d 
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Jet  Concentration  Profiles 


Axial 
locations 
(d  =  ini,  dia.) 

X  =  2d: 


X  =  5d: 


X=  1 0d: 
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AFRL  Water-to- Water  Experimental  setup 


outlet 


Jet  inlet 


Crossflow  inlet 

■=> 


Boundary  conditions: 

Inlets:  Velocities,  species,  K, 
(Fully  devp’d  turb.  profiles) 
Outlet:  Static  Pressure,  K,  z 


Operating  liquid:  Water 


Imaged 

Section 


Laser 

Sheet 


Measurement 


Plane 


Test 

Section 


Collection  _  CCD 
Plenum  *  Camera 


Schematic  of  the  experimental  setup 


Experimental  setup 
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Simulation  matrix 


Considerations: 

•  Turning  of  outlet  Flow 

•  Jets  2,  3,  5,  and  6  not  coplaner  with  exit  pipe 

•  Turbulence  Schmidt  # 

•  Gravity  effects 


outlet 


/v 


Crossflow  inlet 


Plenum  with  exhaust  pipe 


•  Simulation  Matrix: 

•  Single  injector  without  plenum  outlet 

•  Smaller  computational  domain,  less  computational  resources  needed 

•  Single  jet  +  plenum 

•  Multiple  jets:  2,  4,  6 

•  Momentum  flux  ratios:  3.2,  8.2,  14.3,  20.3,  28.3,  36.5,  47.3 
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Grid  Distribution  Around  the  Jet 

(Total  cells:  3  million  Hex  Cells) 
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Grid  Distribution  Around  the  Outlet 
(Total  cells:  3  million  Hex  Cells) 
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^Hybrid  Grid  Around  the  Outlet: 
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Geometry  for  the  Single  jet 
without  the  Outlet  Plenum 


Jet  inlet  (2.36  m/s)  Re  =  11 ,500  (based  on  d) 


i> 


Measure 

< - 

*ment  plane  (x/D  =  3) 

> 

> 

D 

r 

0.63  (m/s)  :  1 .32  (m/s) 
Re  =  25k : 53K 
(based  on  D) 


Grid  at  the  jet  area  (Overall  :2.2  million  hex  cells) 
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Computational  setup  and  Solution  Procedure 


*  Equations  Solved: 

•  Momentum  Eqs. 

•  Continuity 

•  Turbulent  Kinetic  Energy  (TKE) 

•  Turbulent  Dissipation  Rate  (TDR) 

•  Species 

#  Boundary  Conditions: 

•  Inlet  BCs:  Velocities,  K,  z 

•  Inlet  BCs  are  inlet  profiles  obtained  from  separate  fully  developed  pipe 
flow  simulations 

•  Exit  BC:  Static  Pressure 

•  Discretization:  Second  Order  Upwind  for  all  the  equations 

#  Turbulence  Model:  Low  Reynolds  number  K-z  (K-co,  and  other  variants  of  K-z  with 
various  options  were  also  tried) 

•  Material  Type:  Liquid  water 
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Computational  setup  and  Solution  Procedure 


. . .  Continued 

#  Number  of  Grid  Points: 

•  Simulation  performed  on  one  half  of  the  geometry  utilizing  the  gemoetry 
symmetry 

•  Two  grids  generated: 

•  Single  jet  without  plenum:  2.2  million  all  hex  cells 

•  Single  jet  with  plenum:  3  million  mostly  hex  cells 

•  Wall  y+  <  1 

•  Fine  grids  allowed  use  of  the  low  Reynolds  number  K-z  turbulence  model 
integrating  the  solution  to  the  wall 
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Solution  Convergence  History 


Residuals 

continuity 

x-velocity 

-y-velocity 

z-velocity 

k 

epsilon 


—dv- — _ 


Iterations 


Convergence  History  of  Equations'  Residuals 


Convergence  History  of  Mass  Flow  Rate  at  the  Exit 
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Geometry  for  the  Single  jet 
without  Outlet  Plenum 


^Jet  inlet 


Measurement  plane 

Crossflow  inlet 

■=> 

Jet  mass  fraction:  0.05  0.15  0.25  0.35  0.45  0.55  0.65  0.75  0.85  0.95 


Jet  mass  fraction,  J  =  14.3 
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Axial  Velocity  &  TKE  at:  Inlet,  Measurement 
Plane,  and  at  Exit  (J  =  20.3) 


Turbulent  Kinetic  Energy 


Axial  velocity  Turbulent  Kinetic  Energy 


Note  the  Fully  Developed  Turbulent  Inlet  Profile 
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Vortex  Core  and  X-component  of  Vorticity  Distributions 
at  the  Measurement  Plane  and  at  Exit  Plane  (J  =  20.3) 


Jet  inleti 
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Jet  Mass  Fraction  Distributions  &  Penetration 
for  Different  Momentum  Flux  Ratios  (J)  Continued . . . 
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Jet  Mass  Fraction  Distributions  &  Penetration 
for  Different  Momentum  Flux  Ratios  (J) 


J  =  28.3 
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Computed  X-vorticity  &  Jet  Mass  Fraction 
Distributions  at  the  x/D  =  3 
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Computed  and  Measured  Jet  Mass  Fractions 

Continued  .  . 
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Computed  and  Measured  Jet  Mass  Fractions 
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Jet  mass  fraction 


■ 


Computed  and  Measured  Jet  Mass  Fractions 
(Momentum  Flux  Ratio  =  8.2) 


Jet  direction 
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Computational  &  Experimental  Jets 


outlet 


•  Turbulence  Schmidt  number 

•  Computational  domain  does  not  include  the  plenum  and  the  exhaust  pipe 

•  Experimental  Jet  is  not  coplanar  with  exit  pipe 

•  Fluid  issuing  through  experimental  jet  has  to  rotate  &  turn  to  exit 

•  Fluid  issuing  through  computational  jet  has  to  only  turn  to  exit 
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Summary  &  Future  Work 


•  A  series  of  experiments  and  RANS  simulations  are 
performed  to  investigate  JICF  induced-mixing  for  an 
iso-density  liquid  (water). 

•  RANS  simulations,  using  Fluent,  captures  the  flow 
features  such  as  the  CRVs,  and  show  reasonable 
predictions  with  the  measured  concentration  data. 

•  Experimental  data  shows  asymmetry  in  the  wake 

•  Extensive  computations  planned: 

•  Various  momentum  flux  ratios 

•  Multiple  jets  (2,  4,  6) 

•  Plenum  effect 

•  Turbulent  Schmidt  number 
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Computed  and  Measured  Jet  Mass  Fractions 
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Experimental  data  versus  computational  results: 

Axial  Velocity  Contours 


Axial 

locations 

(d  =  inj.  dia.)  Experiment  Fluent  STAR-CCM+  CFD++  LES  (ONERA) 


X  =  5d: 


X=  1 0d: 


II 

i'ii 

III 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Axial 
locations 
(d  =  ini,  dia.) 

X  =  2d: 


X  =  5d: 


X=  lOd: 


LESIONERAi  Ruent 


STAR-CCM+ 


CFD++ 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Experimental  data  versus  computational  results: 


Experiment  Fluent  Star-ccm+  CFD++  LES  (ONERA) 


Axial  Velocity  Contours  and  Radial  Profiles  at  x  «  2d 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


35 


Experimental  data  versus  computational  results: 
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Experimental  data  versus  computational  results: 
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